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The screening length in a quark-gluon plasma, the dispersion relations 
\ of thermal gluon self-energy and the quark potential at high temperature are 

CN 

studied within thermo field dynamics framework. By calculation of the real 
and imaginary parts, of the gluon self-energy in one-loop order in thermo field 

o _ 

\ dynamics, we obtain an expression for the screening length in a quark-gluon 

Qh! plasma and the dispersion relation between the real and imaginary parts. 

D , At high temperature, using photon exchange between electron-positron in a 



X 



skeleton expansion and ladder approximation, the screened Coulomb potential 

is obtained, and using one-gluon and two-gluon exchange between a quark- 

antiquark, we get an expression for the screened quark potential up to 0{g^). 
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I. INTRODUCTION 

The CERN heavy-ion program and the forthcoming start of the BNL-RHIC ( Relativistic 
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Heavy Ion Collider ) make the quark-gluon plasma (QGP) phase of hot QCD of great 
current interest [1]. The screening length and the potential between quarks in a QGP, and 
the dispersion relations of the self-energy of gluon are important topics for the study of 
the QGP phase in a hot QCD. Quantum field theory at finite temperature and density is 
a powerful tool for the study of these subjects. There are three different approaches to 
Quantum field theory at finite temperature and density. It is mainly the following: the 
older Euclidean imaginary-time Matsubara formalism [2] , the closed time-path method due 
to Schwinger and Keldysh [3], and the real time operator-based finite temperature field 
theory, thermo field dynamics (TFD) due to Takahashi and Umezawa [4]. The central idea 
of TFD is the thermal vacuum and the doubhng of the Hilbert space of states. It permits 
thermal averages, which are traces over Fock space, to be written as expectation values with 
respect to the thermal vacuum. The operators in this doubled space are effectively doubled 
in number as well. Within TFD, the quantum field theory is naturally extended to finite 
temperature without having to consider imaginary time. TFD is more concise and an easy 
and convenient way to study the above subjects. 

The studies reported in this paper are based on TFD. Its purpose is three-fold: (1) within 
this framework we calculate perturbatively the screening length in the QGP at one- loop level; 
(2) the dispersion relations of self energy of gluon are derived and, (3) the potential between 
a quark and an antiquark in a hot QCD is obtained up to 0{g'^), g being the QCD coupling 
constant. As a simple example the potential between electron and positron in a hot QED is 
derived. 

The paper is organized as follows: in Sec. II, the screening length in a hot QCD is 
obtained in one-loop level in TFD approach, and we discuss both the color-electric and the 
color-magnetic sectors. In Sec. Ill, the dispersion relations of the gluon self energy at finite 
temperature are proved and explain the physical results on the real part and imaginary part 
of the gluon self energy. In Sec. IV, we calculate the screened Coulomb potential in a hot 
QED up to all orders and the potential between a quark and an antiquark in a hot QCD up 
to 0{g'^). Finally, in Sec. V some conclusions and comments are presented. 
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II. SCREENING LENGTH 



In this section we study the gluon self-energy in one- loop order. To 0{g'^), the gluon 
self-energy is given by the one- loop diagrams in Fig.l. 

The full gluon propagator D'^{x — x') =< A'^{x)Al{x') > may be expressed in terms of 



the one-particle irreducible self-energy, n^5,(/c), 
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where in Minkowski space Q/^i, is the metric tensor with components 

900 = -911 = -922 = -933 = 1; = 0, // 7^ I/. 
Note that D"^^ is diagonal in the color index a and b. U."^^{k) is expressed as 



(1) 



(2) 



(3) 



where U.^^{k) is transverse to one-loop order. 

Before we study the gluon self-energy in one-loop order within TFD framework, it is con- 
venient to rewrite the standard propagators in TFD, A(A;) for bosons and S{p) for fermions, 
as follows [5]: 

for gluons 



A{k)^UBiko,f3)A{k)ui{ko,f3) 



I 



(1 A ^^kij,hir 
(fe2+ie)2 









- 27ri5{k^ 



e/3|fco| _ 1 



£M \ 

1 e 2 



/3|fcol 

e 2 



(4) 



where 



3 



^ 0^ 



(5) 



[/b(A;o,/?) 



^ coshe{ko,(3) smhe{ko,(3) ^ 



smh.9{ko, P) cosh9{ko, P) 



where /3 = ^j^, fc^ and T are the Boltzmann constant and the temperature of the quark- 
gluon system respectively; 
and for quarks 
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where 



S{P) 



%/—m+iST ' 



^ cos^^(b1,/5) e(po)sm^(|p1,/3)^ 
-6(po)sin^(|p1,/?) cos^dpl,/?) 



, sin^e(|p1,/3) = ^, (7) 



and e(po) = 6*^0) - 0{-Pq), 9{po) = 



l,Po > 0, 



0,po < 0. 

The matrices Ub and Up are the Bogohubov transformation matrices [6] for bosons and 
fermions respectively . 
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We calculate the contributions of diagrams in Fig.l to the gluon self-energy in the Feyn- 
man gauge, A = 1, no ghost state appears. In TFD, the operator of the gluon self-energy is a 
(2 X 2) matrix in the Hilbert space of states. Now we express explicitly the (1, 1) component 
of the contributions of the diagram 1(a) to the gluon self-energy. 



^ f J (27r)4^^^'^ [{p - kY -w? + ie] (p2 - m2 + ie) ^ 



{2tt 

where Tr{T"-T^) — QS"^, T", are generators of the non-abelian gauge group, Nf is the 
quark flavor number, and Q is the Casimir operator related to the fermion representation of 
the non-abehan gauge group. We normahze the generators T" and T'^, so that Q — 1. The 
first term on the right side of Eq. (8) is independent of temperature, T — 0. The second 
and third terms in the right side of Eq. (8) depend on temperature, T ^ 0. The second 
term is the real part of the gluon self-energy and the third term is the imaginary part of 
the gluon self-energy. In order to obtain the pole mass of a gluon at finite temperature, we 
consider the (0, 0) component of 11^^ and let the momentum /c — > and the energy ko — >■ 0. 
For simplicity, we omit the quark mass m in Eq. (8). Then, we get the contribution of the 
diagram (a) in Fig.l to the square of the pole mass of a gluon at finite temperature from 
Eq. (8) 

Relllt(k ^ 0, fc„ ^ 0, r ^ 0)];,'.„„., = /~ (9) 

where in the integral for pq, we have used this formula 
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5(p-^-m') = —[6{po + u;) + S{po-u;)], uj = ^p' + m\ (10) 



After using the following integral 

POO ^cjT 1 

/ ;^^ = 7;^(i-2'"'^)r(^K('^)' (i?e/. > o, i?e^ > o) (ii) 

where and C,{v) are the Gamma function and the Riemann's zeta function respectively, 

we get 

i?e[nSS(^^ 0,fco - 0,T ^ 0)]i.\,.,(„) = ^^'^"^' ^^r(2)C(2) (12) 

Similarly we have obtained contributions of diagrams (b) and (d) in Fig.l to the square 
of the pole mass of a gluon at finite temperature, as follows: 

Re[nf,{k ^ 0, A;o ^ 0, T ^ ^)]%mM) = ^^^r(2)C(2) (13) 

where 

E Co^cCm = C5^^ C^bc = -Cbac; C = TV, /or 5[/(7V) ^aii^e ^roi^p. (14) 

a,b 

Here Co^c are structure constants of the non-abelian gauge group, and C is the Casimir 
operator related to the adjoint representation of the non-abelian gauge group. 
In the calculation of Eq. (13) we have used the following formula: 

r^^^;?^^^^)^^'^) (i?e/.>0,i?ei.>l). (15) 

To sum up, we obtain finally the total contributions of all diagrams in Fig.l to the square 
of the pole mass of a gluon at finite temperature in one-loop order 



m. 



^-^(iv + -iv,)r(2)c(2) 

S^''l{N + ^)g'T^ (16) 



where 



OO -| 

r(2) = 1, c(2) = E - 




(17) 



and the mass rriei is the color electric mass of a gluon at finite temperature in one-loop order. 

By calculations for the spatial components of the gluon self-energy in Fig.l in TFD, it 
is found that U^^{ko,k ^ 0,T ^ 0),n«^ and T,f=i'^u are all zero. So, the color magnetic 
mass of a gluon at finite temperature in one- loop order equals zero. Therefore in one- loop 
order the color-electric gluon fields are screened and the color-magnetic gluon fields are not 
screened. The conclusions from the present calculation of the gluon self-energy in one-loop 
order in TFD are the same as in the conventional imaginary-time formalism (ITF) [7] , but 
the calculations in TFD are more clear, simple and easy. The screening length in a QGP is 
the reciprocal of niei, then 



In this section we discuss what is the dispersion relation for the thermal gluon self- 
energy? For this purpose, we consider the renormalized gluon thermal propagator 



a/3 



(18) 




III. DISPERSION RELATION 



(19) 



where the Feynman gauge has been chosen, A = 1, and n(A;) and Il{k) are connected with 




each other through the following relation 



W'^ik) = [TUB{ko)n{k)UB{ko)T]''^ 
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(21) 



where the i^) vector indices are omitted and 



n"/3 ^ r'^'^ReU + iS'^'^ImU. (22) 
Exphcitly, the (1,1) component of the relation (20) is 

11 _ cosh^^(/co) smh^ejko) 

[^R{k)U - 3.-^\ko + irr-u;^ ~ {ko-iVr-J ^^^^ 

uj^^p + rn^ + ReU - r^, ImU = -2kor (24) 

From Eqs. (22) and (24) it is clear that the real and imaginary part of our interest are given 
by it!en and ImU. Prom Eqs. (21) and (22) we obtain 



ReU^^ = -i^e^2^ ReU^'^ = ReU^^ = 0, (25) 



ImW^ = ImW, ImU^' = ImU'\ (26) 



1 pf^\ko\ I 1 

^-n- = -i^-^M-. (27) 



2 e/^l'^ol/^ 

Further, we get 



ReU = i?en^\ (28) 

p/3|fco| _ 1 

ImU = ImU^^ (29) 

e'^lfcol 

/mn = -rr^ImU^'^. (30) 

These equations for i^ell^^ and Jmll^^ are interesting. It is clear that U^^ and 11^^ arc purely 
imaginary, and the real part of the gluon self-energy can be obtained from U^^ or H^^ [8]. 
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By calculating U.^^{k), for ^ ^ 0, and keeping the energy kg ^ 0, we obtain the imagi- 
nary part of the (0, 0) component of the gluon self-energy for the diagram (a) in Fig.l 

7mn(A; = 0, fco) = -^^-^^lm^'\k^) 



27r 2e'3l'=ol/2 {emM-^ + 1)^ 
g^QS'^^Nf (e'^l^^ol/^ - l)kl 



(31) 



47r (e/3|feo|/2 + 1) 
where we have omitted the (0, 0) vector index. 

By calculating 11^^ (A;), for ^ — > 0, and keeping the energy k^ ^ 0, we obtain the real part 
of the (0, 0) component in the vector indices of the gluon self-energy for the diagram (a) in 
Fig.l 



ReYi{k = Q,ko) = ReW\k 



0) 



(32) 



7r2 Jo (e/5- + l)(A;g-4cu2) 
where we omitted also the (0, 0) vector index. From Eqs. (31) and (32), we get 



RelLiko) = - r ^^^^sinh^(^X (33) 

TT J-oo kn — kn 2 



-OO AtQ — AtQ 

This is the dispersion relation of the gluon self-energy in Fig. 1(a). 



IV. QUARK POTENTIAL 

In this section we study the potential between a static quark and anti-quark within TFD. 

We consider the free energy of a configuration of Ng quarks and Ng anti-quarks, 
F(fi, f2, . . . rW^; ri, r2 • • • ^n^)- % definition 



9 



where H is Hamiltonian of the system and in the state, \S), there are Ng static heavy quarks 
at fi, ■■■f]^q and Nq static heavy anti-quarks at r^, ■■■rj^^, and Z]|s)(5'|5') = N. 
We can write 



Y^iSle-^^'lS) = Y.{S'\ E ^a,(ri)V'a.(rl)...V'a.,(rlvjV'6^,(r1)V',^,(r",)...V',V(4Je-'^'' 

\s) W) {a,b) 



< (rM, ir2)..4l^ {tnMI (rMl (^1)-<- (4,) I '5') (35) 

where these quarks and anti-quarks are taken from \S), and in \S') there are no quarks 
and anti-quarks, and V'ai(ri) ('0li(^i)) ^-iid V'ftil^i) ('06i(^i)) ^'I'e annihilation ( creation ) 
operators of quark and anti-quark respectively. 
Prom Eqs. (34) and (35) we get 

e-PF^,,N, ^ Tr[e-'^^L(fi)...L(r-WjLt(r-;)...Lt(4^_)] (36) 

where the trace, Tr, is J2\s')y only over states of the pure gauge field and includes no quark 
field. The quarks are considered as external particles, and the Wilson loop, L{r), is 

L(f) = lTr[Te/o^^"^S(-»'='-] (37) 
where T represents time ordering. 

For the quark operator, ■0, we have 

if-P-^l- (38) 

H — —igi/j^ij,Tt/jAn + rmjjip — —igijj^oTiljAo + rm/ji/j (39) 
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due to the quark being static, with the current, J — and Jq ^ 0. 
Using the commutation relations 



bPiirl, t),Mrj, t)]+ = SiAb, (40) 

we have 

[H, -0]- = -igrAoip + mip. (41) 

We put m — 0, hence we get 



]^^r-A-M^{x). (42) 

Prom Eq. (42), we get 



Mr, t)^T exp[ f gT'^Alir, T)dT^a{r, 0) (43) 

Then the relation given in Eq. (36) can be derived from Eq. (43). For one static quark, Eq. 
(36) becomes 



g-/3Fi,o ^ Tre-'^^L(f) (44) 



and 



g-/3Fo,o ^ Tre"/^^. (45) 



Hence we have 
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where Fo,o is the free energy of the pure gluon field and Fi_o is the free energy of one quark 
and the gluon field. Since 



Tr[e-'^^L(f)] = ^(s'|e-^^L(f)|s'), (47) 

\s') 



Tre-'^^ = X^(s'|e-'5^|s') (48) 



where the \s') is the state with only gluons and 



H\s') = Es>\s'), (49) 



we have 



We can scale Egi up a times and continue to the complex plane, and given that 



(50) 



a = pe'\ e = -e, (51) 
and put Rea — > oo and Ima < 0, then —a/3 — > — oo, therefore we get 



e"^^^ = (L(r)) (52) 

where Fg is the free energy of a single-quark system without the quanta of pure gauge field. 
For a static quark and anti-quark we have 



-l3V{r-r') 



= {L^{f)L{P)), 
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(53) 



where V{f— r) is the potential between a static quark and anti-quark, and the lower index 
c representing the connected two-point function of the Wilson loop. 

First we consider the case of the abelian gauge field. This is easily done in QED, by 
including the contributions of all Feynman diagrams of virtual photons exchanged between 
an electron and a positron. Now we consider the ee potential energy at finite temperature. 
In this case, starting with a skeleton expansion, in which only full photon propagator and 
the renormalized coupling constant e are used, then the screening potential is the summation 
of a set of ladder graphs shown in Fig. 2. 
Therefore we get 

g-/3VM) = 1 + y lr^g-m(/3)r]n _ lll^-m{p)r^ ^r,^-. 

^^n\ATrr Airr 
This simply leads to a screened Coulomb potential 

V(r,p) = -—e-'^^^^, (55) 
47rr 

where m(/3) is the renormalized screening mass of a photon at finite temperature. The 
magnitude of m(/3) is given by an expression similar to Eq. (16). 

Now we study the case of QCD in TFD. For the connected two-point function of Wilson 
loop in Eq. (53), up to 0((?^) the contributions are from the diagrams in Fig.3. 

The contribution of the one-gluon exchanged diagram 3(a) to the connected two-point 
function of Wilson loop equals zero since 



Tr(T") = 0. (56) 

For the two diagrams 3(6) and 3(c), using the thermal gluon propagator in Eqs. (4) and 
(5), we have obtained the contribution of 0{g^) to the connected two-point function of the 
thermal Wilson loop as follows: 
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where Tr(T'^T^) = Q5°-^ , g is the coupling constant, T is the temperature, and the factor 
(A^^ — 1) is from the adjoint representation of the gluon for the gauge group SU{N). Here 
we have used the following formulas: 



Jo (l + e^)2 ^i(a + A;)"' ^ ^ 

which has the value 



k=l 

for a = and n = 2 
and 



where S(| — /3, | — a) is the Beta function. 

The function G(p), in Eq. (57), can be Fourier transformed to the coordinate space, as 
follows: 



2aU^ (61) 



If we take into account the thermal gluons in Fig. 3 with the color-electric pole mass, 
meiiT) and use the following formula 
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/ (k: 



An' 



-2m^l{T)\x-x I 



(F + ml{T))[{p- ky + ml{T)\ {x - xj 



(62) 



we get 



g^Q^{N^ - 1) Q-^rn,m\s-s'\ 



G{x-x)^ ^^^^^ ^,.2 • (63) 



Therefore the potential between a quark and an anti-quark at high temperature up to 0{g^) 
is 



V(r) = ^ ^ ^ . (64) 

It is clear from Eqs. (55) and (64) that the potential in a QCD plasma is different than 
the potential in a QED plasma, though they all decrease exponentially with the increase of 
the distance, r, between them, but the decrease in power are different, the latter is linear 
while the former is quadratic [9]. 



V. CONCLUSIONS AND DISCUSSIOS 

In summary, we have obtained several consequences for QGP within TFD framework: 
the screening length in QGP in one-loop order, the dispersion relations between the real 
part and the imaginary part of the thermal gluon self-energy and the quark potential at 
high temperature up to 0{g'^). The main advantage of TFD compared to other formalisms 
of thermal field theories is the ease of calculation. 

Hopefully QGP will be observed by experiments at CERN and at BNL-RHIC in the 
near future. As an estimate, the binding radius Rb of J/ijj particle [10] can be considered 
to be roughly Rb ~ 0.3/m [11], setting the running coupling constant g{T) ~ 1.0, then, we 
obtain the solution L(T) — Rb and the critical temperature Tc ^ O.SGev. This means that 
a suppression of J/ip mesons should be observed if the temperature in a nuclear collision 
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becomes much higher than the transition temperature Tc- 
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Figure 1: One-loop self-energy diagrams of a gluon. The wavy line represents 
the gluon, the solid line represents the quark or anti-quark and the dotted 
line represents the ghost. The loop diagram with shaded area represents the 
one-loop self-energy of a gluon. 




Figure 2: The ladder graphs of photons exchanged between an electron and 
a positron. 




(a) (b) (c) 



Figure 3: The connected two-point function of the thermal Wilson loop up 
to O(g^). The wavy line represents the internal thermal gluon, and the 
quark and anti-quark are considered as external particles. 
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